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A Simple Explanation for the X(3872) Mass Shift Observed for Decay to D*^D^ 
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We propose a simple explanation for the increase of approximately 3 MeV/c^ in the mass value of 
the X{3872) obtained from D'^D'^ decay relative to that obtained from decay to J/ip-K^-K~ . If the 
total width of the X(3872) is 2-3 MeV, the peak position in the D'^D'^ invariant mass distribution 
is sensitive to the final state orbital angular momentum because of the proximity of the X(3872) to 
D*Oj=jo threshold. We show that for total width 3 MeV and one unit of orbital angular momentum, a 
mass shift ~3 MeV/c^ is obtained; experimental mass resolution should slightly increase this value. 
A consequence is that spin-parity 2^ is favored for the X(3872). 

PACS numbers: 14.40.Gx, 13.25.Hw 



The X(3872) has been seen primarily in its J/-07r+7r 
decay mode P, 0, S, 0, S], from which a measured 
mass value of 3871.4 ± 0.6 MeV/c^ is obtained 0; an 
upper limit on the total width of 2-3 MeV has been 
estimated [l[. Observation of decay to the J/'07 fi- 
nal state U Q has established positive C-parity, and 
analysis of the decay angular distributions |9| has nar- 
rowed the spin-parity (J^) possibilities to 1"*" or 2~ for 
the X(3872). The invariant mass distribution for the 
D^D^tt'^ system resulting from B-decay to the D^D'^tt^K 
final state shows a peak near threshold yielding a mass 
value 3875.2 ± 0.7 + 0.9 - 1.8 McV/c^ [l^l, and this 
has been interpreted as evidence for the decay process 
X(3872) D*°D°. A recent BaBar analysis [lH has 
confirmed this result, and has obtained the correspond- 
ing mass value 3875. I^q'^ ± 0.5 MeV/c^. In each case, 
the first error quoted is statistical, and the second is 
systematic. These results are in excellent agreement. 
Compared to the mass value 3871.4 ± 0.6 MeV/c^ 0] for 
the J/i/'TT+TT^ decay mode, the mass difference is 3.8^2 q 
MeV/c2 from Belle, and 3.7to,l MeV/c^ from BaBar. 
The significance of the latter is at the four standard devi- 
ation level, and given the consistency of the BELLE and 
BaBar results for D*^D''\ it would seem to be a real effect. 
We take the point of view that this is indeed the case, 
and suggest a possible explanation, which, although very 
simple, carries some significant physical implications. 

Consider the decay process B KX, with X is a 
resonance decaying to a final state F. The invariant mass 
distribution for the system F takes the form 
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where Ci is a constant, m is the invariant mass of system 
F, TB(rn) is the invariant amplitude describing the B to 
KX coupling, Tpim) is the invariant amplitude describ- 
ing the X to F coupling, d(j)F{m) is the element of F de- 
cay phase space, and the denominator is the square of the 
relativistic Breit-Wigner propagator describing the reso- 
nance X. The factor m is present because the Lorentz- 
invariant phase space volume element is proportional to 



dm^, and (p/ms) is the two-body phase space factor for 
B KX decay, with 



Equation (1) can be written in terms of the partial 
width for X decay to F, rp{m), as 



dN ^ C2-m-p- |rB(m)prf(m) 
dm (ml-m^f +mlrtot{my' 

In general, Ttot{fn), which is the mass-dependent total 
width of X, takes the form 
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rtot(m) = y^ri(m), 



(4) 



where M is the number of decay modes of X, and the Ti 
are the individual partial widths, of which F^? is one. 

The amplitude Tb is not known. However, if X has 
spin J, angular momentum conservation requires that 
there be J units of orbital angular momentum associ- 
ated with the KX system resulting from B decay, and so 
we express Tb in terms of the corresponding centrifugal 
barrier factor as follows: 



TB{m) 



(5) 



where Dj{p,R) is the Blatt-Weisskopf Damping Fac- 
tor [l^ . and R is the associated radius parameter, for 
which we choose the value 5 GeV~^ (i.e. 1 Fermi). The 
Dj functions for J = — 3 are summarized in Table 
I. For the X(3872), the mass range is limited (3.87-3.91 
GeV/c^) and the Q-valuc for B KX decay is large 
( 0.9 GeV/c^), so that the m-dependence introduced by 
this description of the B decay vertex is small. 

For the X(3872ji onlythree decay modes have been ob- 
served to date [H, B, B 11, S S S, Efl, El , and of these the 



2 



J Dj{p,R) 



1 

1 1 + {pRf 

2 9 + 3(pi?)' + (pi?)* 

3 225 + 45(pfl)^ + &{pRY + (pi?)*^ 

TABLE I: The Blatt-Weisskopf Damping Factors. 

J/ipTT~^TT~ and J/tpl modes are well above their respec- 
tive mass thresholds. It follows that the m-dependence 
for each of these contributions to Eq.(4) is very weak. 
In contrast, for the D*°D'^ mode, the dependence on m 
is potentially quite strong because of the proximity of 
the invariant mass threshold to the X{3872) mass; this 
is, in fact, the main point of this Letter. If the state F 
of Eq.(3) is chosen to be the D*^D^ final state, the nu- 
merator is very sensitive to its dependence on m, as will 
be discussed in detail below. However, in the denom- 
inator, this mass dependence is subsumed into that of 
the total width by way of Eq.(4). Since the other known 
modes have very little mass-dependence, and since other 
as yet unobserved decay modes may contribute to Ttot, 
we choose to treat r^ot a-s a constant in the computations 
to follow. 

From Eqs.(l) and (3), 
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TF{m) ^ / |Tf(to)| d0F(m), 



and with F the D*'^D" final state. 



(7) 



where 



2^ 

is the momentum in the D*^D^ rest frame; we ignore the 
width of the D*^ (< 1 MeV [1]). Integrating over f], we 
obtain 



TFim) {q/m)\TF{m)\ 



(9) 



and following Eq.(5) (with the same value of R) we ex- 
press Tp{ra) as 



VDL{q,R) 



(10) 



where L is the number of units of orbital angular momen- 
tum associated with the D*^D^ system resulting from 
decay of the X (3872). 

The full expression for the D*^D^ invariant mass pro- 
jection is then 

diV ^ {p^''+^/Dj{p,R)) {q'^+'/DL{q,R)) 
-; — — : ^ : — : (ii) 



dm 



(m^ - 777,2)^ +TO§r2„t 



where C3 is a constant. 

If the X(3872) has ^ 1+ L can be or 2, while for 
= 2^ L can be 1 or 3, since parity should be conserved 
in the decay process. Consequently we consider that only 
L values in the range 0-3 are of relevance to the X(3872). 
Also we use the following central mass values @ in our 
calculations: 

m{D°) ^ 1864.84 [±0.18] MeV/c^ 
m{D*°) = 2006.96 [±0.19] MeV/c^ 



These yield the threshold_ mass value 3871.80 ± 0.37 
McV/c^ for decay to D*^D'^ (again ignoring the width 
of the D**^), where the error is obtained by combining 
the error on twice the mass and that on the D*^ — 
mass difference in quadrature. 

We choose X(3872) mass values 3870.8, 3871.4 and 
3872.0 MeV/c^ (i.e. the PDG 2007 average and plus 
or minus one sigma), and for each use Ttot values 2, 3 
and 4 MeV. For each of these nine combinations we use 
Eq.(ll) to compute the lincshapc for the choices L ~ 
and 2 with J = 1 (i.e. X(3872) = 1+), and L = 1 and 
3 with J = 2 (i.e. A:(3872) = 2"). For L = 2, all of 
the lineshapes obtained are very broad, reaching a maxi- 
mum close to mass 3.90 GeV/c^, while for L = 3 the dis- 
tribution obtained increases monotonically through 3.91 
GeV/c^, the upper limit of the region investigated. It 
follows that for L = 2 and L = 3 Eq.(ll) yields behavior 
which is totally unlike that observed for data [13, [lH, 
and so the possibility that such contributions play a sig- 
nificant role in X(3872) decay to D*^D^ is discarded. 

For the remaining [J^ = 1+ : L = 0, J = 1) and 
(J^ = 2~ : L = 1, J = 2) possibilities, Eq.(ll) does lead 
to lineshapes which peak at mass values a few MeV/c^ 
larger than the input X(3872) values. The mass shift is 
defined as the difference between the observed peak mass 
position and the input X(3872) mass value, mo, and the 
mass shifts obtained are summarized in Tables II and III 
respectively. In both Tables, the mass shift increases with 
decreasing mo and increasing Ttot- However each value 
in Table II is smaller than the corresponding one in Table 
III, the difference averaged over the Tables being ~1.85 
MeV/c^. Furthermore, it seems from Table II to be un- 
likely that a mass shift of ~3 MeV/c^ can be obtained 
for a reasonable choice of AT (3872) mass and width, es- 
pecially since the D*'^D^ mass distribution which results 
becomes considerably broader as the X(3872) mass is 
decreased and its width increased. 

Using X(3872) mass 3871.4 MeV/c^ and width 3 
MeV, we illustrate the lineshape behavior obtained from 
Eq.(ll) for JP = 1+ and .F = 2" in Fig.l and Fig.2 
respectively. In each figure, the curve is obtained us- 
ing Eq.(ll), and is then used to generate the 3000 
events shown in the histogram below, which uses the 
mass intervals from the BaBar analysis [llj. The peak 
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X{3872) mass [MeV/c^] 


2 


3 


4 


3870.8 


1.54 


1.75 


1.99 


3871.4 


0.90 


1.17 


1.45 


3872.0 


0.45 


0.75 


1.02 



[MeV] 



TABLE II; Dependence of the peak mass shift (in MeV/c^) 
on X(3872) mass and total width for L = and J — 1 (i.e. 
= 1+). 



X(3872) mass [MeV/c^] 


2 


3 


4 


3870.8 


3.79 


4.25 


4.75 


3871.4 


2.37 


2.99 


3.61 


3872.0 


1.28 


1.98 


2.66 



Ftot [McV] 



TABLE III: Dependence of the peak mass shift (in MeV/c^) 
on X(3872) mass and total width for L = 1 and J = 2 (i.e. 



shift in Fig. 2 (a) agrees well with the result from experi- 
ment 11 1 , and the observed signal shapes seem better 
represented by that of Fig. 2(b) than by that of Fig. 1(b). 
However, it must be acknowledged that the experimen- 
tal uncertainties are significant, and that even the uncer- 
tainty in the location of the D*^D^ threshold (above) is 
relevant on the scale of the effect under discussion. 
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FIG. 1: (a) The m[D*°D°) lineshape obtained from Eq.(ll) 
for J-^ = 1+ using X(3872) mass 3871.4 MeV/c^ (indicated 
by the dot-dashed line) and width 3 MeV. (b) The histogram 
for 3000 events generated using the curve shown in (a). 
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FIG. 2: (a) The m{D*°D") lineshape obtained from Eq.(ll) 
for = 2" using X(3872) mass 3871.4 MeV/c^ (indicated 
by the dot-dashed line) and width 3 MeV. (b) The histogram 
for 3000 events generated using the curve shown in (a). 



Near threshold, such effects should not be represented 
by Gaussian smearing in mass, as is done usually, since 
this will yield contributions below threshold. It is three- 
momentum resolution which is the source of the smear- 
ing, and this must be investigated by full detector sim- 
ulation for the experiment in question. Since such sim- 
ulation obviously cannot yield events below threshold, 
it seems probable that the peak mass shifts calculated 
above will be increased as a result of experimental res- 
olution. We suspect that such effects will be small (< 
1 MeV/c^), but a thorough investigation making use of 
detector simulation is necessary. 

In summary, we have shown that a simple treatment 
of the orbital angular momentum involved in X(3872) 
decay to D*^D^ can account for the difference between 
the mass measured in this mode and that obtained from 
J /ipTT^n' decay. The results favor ~ 2^ over = 
l"*" for the X(3872), but the uncertainty in the measured 
mass difference (3.7 ±1.2 MeV/c^), and the absence of 
simulated detector resolution effects, prevent a definite 
conclusion. If our interpretation is correct, a corollary is 
that the width of the X(3872) cannot be much smaller 
than ^2 MeV, since otherwise significant displacement 
of the invariant mass peak for D*'^D^ would not occur. 

Work supported by the U.S. Department of Energy 
under contract number DE-AC03-76SF00515. 



We have made no attempt to study the effect of de- 
tector resolution on the mass shifts calculated above. 
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